arXiv:1505.01223vl [cond-mat.supr-con] 5 May 2015 


Direct evidence for the emergence of a pressure induced nodal 
superconducting gap in the iron-based superconductor 


Bao.65R't>o.35Fe2As2 


Z. Guguchia/’F] A. Amato , 1 J. Kang , 2 H. Luetkens , 1 P.K. Biswas , 1 
G. Prando , 3 F. von Rohr , 4 Z. Bnkowski , 5 A. Shengelaya , 6 
H. Keller , 4 E. Morenzoni , 1 R.M. Fernandes , 2 and R. Khasanov 1 

laboratory for Muon Spin Spectroscopy, 

Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland 
2 School of Physics & Astronomy, University of Minnesota, Minneapolis, MN 55455, USA 
3 Leibniz-Institut fur Festkorper- und Werkstoffforschung 
(IFW) Dresden, D-01171 Dresden, Germany 
Physik-Institut der Universitat Zurich, 

Winterthurerstrasse 190, CH-8057 Zurich, Switzerland 
5 Institute of Low Temperature and Structure Research, 

Polish Academy of Sciences, 50-f22 Wroclaw, Poland 
6 Department of Physics, Tbilisi State University, 

Chavchavadze 3, GE-0128 Tbilisi, Georgia 


1 


Abstract 

Identifying the superconducting (SC) gap structure of the iron-based high-temperature super¬ 
conductors (Fe-HTS’s) remains a key issue for the understanding of superconductivity in these 
materials. In contrast to other unconventional superconductors, in the Fe-HTS’s both d-wave and 
extended s-wave pairing symmetries are close in energy, with the latter believed to be generally fa¬ 
vored over the former. Probing the proximity between these very different SC states and identifying 
experimental parameters that can tune them, are of central interest. Here we report high-pressure 
muon spin rotation experiments on the temperature-dependent magnetic penetration depth A (T) 
in the optimally doped Fe-HTS Bao.65Rbo.35Fe2As2. At ambient pressure this material is known to 
be a nodeless s-wave superconductor. Upon pressure a strong decrease of A (0) is observed, while 
the SC transition temperature remains nearly constant. More importantly, the low-temperature 
behavior of 1/A 2 (T) changes from exponential saturation at zero pressure to a power-law with 
increasing pressure, providing unambiguous evidence that hydrostatic pressure promotes nodal SC 
gaps. Comparison to microscopic models favors a d -wave over a nodal s" 1 -wave pairing as the 
origin of the nodes. Our results provide a new route of understanding the complex topology of the 
SC gap in Fe-HTS’s. 

PACS numbers: 
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I. INTRODUCTION 


After six years of intensive research on the Fe-based high temperature superconductors 
(Fe-HTS’s), no consensus on a universal gap structure has been reached. There is evidence 
that small differences in electronic or structural properties can lead to a strong diversity in 
the superconducting (SC) gap structure. On the one hand, nodeless isotropic gap functions 
were observed in optimally doped Bai_ a ,K T Fe2As2, Bai_ x Rb 3 ; Fe2As2 and BaFe2- x Ni x As2 as 
well as in BaFe2-xCo x As 2 , K a .Fe 2 _ ?; Se2, and FeTei_a,Se x nm. On the other hand, sig¬ 
natures of nodal SC gaps were reported in LaOFeP, LiFeP, KFe 2 As 2 , BaFe 2 (Asi_, E P ;E )2, 
BaFe2_a;Rux-As2, FeSe as well as in over-doped Bai_ x K r Fe2As2 and BaFe2- a; Ni ;c As2 [ 3 E 0 - 
□3J. Understanding what parameters of the systems control the different SC gap structures 
observed experimentally is paramount to elucidate the microscopic pairing mechanism in 
the Fe-HTS’s and, more generally, to provide a deeper understanding of the phenomenon 
of high-temperature superconductivity. On the theoretical front, it has been proposed that 
both the s +_ -wave and d -wave states are close competitors for the SC ground state [T%H 2 Tj . 
Although the former generally wins, it has been pointed out that a d -wave state may be 
realized upon removing electron or hole pockets. On the experimental front, a sub-leading 
d -wave collective mode was observed by Raman experiments inside the fully gapped SC 
state of optimally doped Bai_ a ,K r Fe2As2 [251 [ 26 ] . In KFe2As2, a change of the SC pairing 
symmetry by hydrostatic pressure has been recently proposed, based on the U-shaped pres¬ 
sure dependence of T c [ 27 ]. However, no direct experimental evidence for a pressure induced 
change of either the SC gap symmetry or the SC gap structure in the Fe-HTS’s has been 
reported until now. Here, we show unambiguous evidence for the appearance of SC nodes 
in optimally-doped Bai_. r Rb x Fe2As2 upon applied pressure, consistent with a change from 
a nodeless s +_ -wave state to either a nodal s +_ -wave or a d -wave state. 

Our results rely on measurements of the magnetic penetration depth A, which is one 
of the fundamental parameters of a superconductor, since it is related to the superfluid 
density n s via 1 /A 2 = //o e 2 n s /m* (where m* is the effective mass). Most importantly, the 
temperature dependence of A is particularly sensitive to the presence of SC nodes: while 
in a fully gapped SC AA -2 (T) = A ~ 2 ( 0 ) — A -2 (T) vanishes exponentially at low T, in a 
nodal SC it vanishes as a power of T. The muon-spin rotation (//SR) technique provides a 
powerful tool to measure A in type II superconductors [ 2 E] • A //SR experiment in the vortex 
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state of a type II superconductor allows the determination of A in the bulk of the sample, 
in contrast to many techniques that probe A only near the surface. 

For the compound Ba a 65 Rbo.35Fe 2 As2 investigated here, and for the closely related system 
Bai_ x K 2 .Fe2As2, previous pSR measurements of A (T) revealed a nodeless multi-gap SC state 
M, in agreement with ARPES measurements p Eli m\. In this communication, we 
report on pSR studies of A ( 0 ) and of the temperature dependence of AA ~ 2 in optimally 
doped Bao.65Rbo.35Fe 2 As 2 under hydrostatic pressures. This system exhibits the highest T c 
~ 37 K among the extensively studied “ 122 ” family of Fe-HTS’s. We observe that while T c 
stays nearly constant upon application of pressure, A ( 0 ) decreases substantially. In view of 
previous works in another “ 122 ” compound that reported a sharp peak of A( 0 ) at a quantum 
critical point EH. we interpret the observed suppression of A( 0 ) as evidence that pressure 
moves the system away from a putative quantum critical point in Bao.65Rbo.35Fe2As2. More 
importantly, we find a qualitative change in the low-temperature behavior of AA -2 (T) 
as pressure is increased. While at p = 0 an exponential suppression characteristic of a 
nodeless superconductivity is observed, for p = 2.25 GPa a clear power-law behavior is 
found. Because pressure does not affect the impurity concentration, which could promote 
power-law behavior even for a nodeless system [3 23, our findings provide strong evidence 
for a nodeless to nodal SC transition. Our fittings to microscopic models reveal that this 
behavior is more compatible with a d -wave state rather than an s +_ state with accidental 
nodes, indicating that pressure promotes a change in the pairing symmetry. 


II. RESULTS 

A. Probing the nonuniform field distribution in the vortex state under pressure 

Figures la and b exhibit the transverse-field (TF) /iSR-time spectra for 
Ba 0 .65Rbo.35Fe 2 As2, measured at ambient p = 0 GPa and maximum applied pressure p 
= 2.25 GPa, respectively. The spectra above (45 K) and below ( 1.7 K) the SC transition 
temperature T c are shown. Above T c the oscillations show a small relaxation due to the ran¬ 
dom local fields from the nuclear magnetic moments. Below T c the relaxation rate strongly 
increases with decreasing temperature due to the presence of a nonuniform local magnetic 
field distribution as a result of the formation of a flux-line lattice (FLL) in the SC state. 
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Figure 1: (Color online) Transverse-field (TF) ^SR time spectra and the corresponding 
Fourier transforms (FT’s) of Bao.65Rbo.3sFe2As2. The spectra are obtained above (45 K) 
and below (5 K) T c (after field cooling the sample from above T c ): (a,c) p = 0 GPa and (b,d) p 
= 2.22 GPa. The solid lines in panels a and b represent fits to the data by means of Eq. 3. The 
solid lines in panels c and d are the FT’s of the fitted time spectra. The dashed and solid arrows 
indicate the first moments for the signals of the pressure cell and the sample, respectively. 

Figures lc and d show the Fourier transforms (FT’s) of the /rSR time spectra shown in 
Figs, la and b, respectively. At T = 5 K the narrow signal around poH ext = 50 mT (see 
Figs, lc and d) originates from the pressure cell, while the broad signal with a first moment 
/loFfint < PoHext, marked by the solid arrow in Fig. lc, arises from the SC sample. 

Below T c a large diamagnetic shift of HoH int experienced by the muons is observed at all 
applied pressures. This is evident in Fig. 2a, where we plot the temperature dependence of 
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Figure 2: (Color online) Diamagnetic shift A£>dia (a) and the muon spin relaxation rate 
<7 SC (b) of Bao. 65 R.bo. 35 Fe 2 As 2 as a function of temperature at various pressures, (a) The 

definition of the diamagnetic shift Ah?di a is given in the text, (b) The muon spin relaxation rate 
<t sc is measured in a magnetic field of poH = 50 mT. The dashed vertical lines denote T c for p = 
0 and 2.22 GPa. 



the diamagnetic shift AF>di a = ho[-^int,sc--^mt,Ns] for Bao. 65 Rb>o. 35 Fe 2 As 2 at various pressures, 
where ^o-^int,sc denotes the internal field measured in the SC state and /io^int,NS the internal 
field measured in the normal state at 45 K. Note, that poH int ^s is temperature independent. 
This diamagnetic shift indicates the bulk character of superconductivity and excludes the 
possibility of field induced magnetism [33j in Bao. 65 R.bo. 3 sFe 2 As 2 at all applied pressures. The 
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Figure 3: (Color online) Pressure evolution of A ~ 2 (T) of Bao. 65 Rbo. 35 Fe 2 As 2 . The tempera¬ 
ture dependence of A -2 measured at various applied hydrostatic pressures for Bao. 65 Rbo. 3 sFe 2 As 2 . 
The solid line for p = 0 GPa corresponds to a two-gap s-wave model (a) and the solid lines for finite 
pressure represent a fits to the data using a multiband d -wave model (b-f). The dased lines in pan¬ 
els (c) and (f) represent fits to the data using the microscopic model described in the supplemental 
material. 

SC transition temperature T c is determined from the intercept of the linearly extrapolated 
Ai?dia curve its zero line (we used the same criterium for determination of T c from ABfc a (T) 
as from the susceptibility data / \ m (T), presented in the supplementary material). It is found 
to be T c = 36.9(7) K and 35.9(5) K for p = 0 GPa and 2.25 GPa, respectively. The ambient 
pressure value of T c is in perfect agreement with T c = 36.8(5) K obtained from susceptibility 
and specific heat measurements (see supplemental material). At the highest pressure of p — 
2.25 GPa applied, T c decreases only by ~ 1 K, indicating only a small pressure effect on T c 
in Ba 0 . 65 Rbo. 35 Fe 2 As 2 . The temperature dependence of the muon spin depolarization rate 
cr sc , which is proportional to the second moment of the field distribution (the moments of 
the field distribution probed by the muons were extracted with the equations described in 
the Method section), of Bao. 65 Rb 0 . 35 Fe 2 As 2 in the SC state at selected pressures is shown 
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in Fig. 2b. Below T c the relaxation rate a sc starts to increase from zero with decreasing 
temperature due to the formation of the FLL. It is interesting that the low-temperature 
value <t sc (5 K) increases substantially under pressure (see Fig. 2b): cr sc (5 K) increases about 
30 % from p = 0 GPa to p = 2.25 GPa. Interestingly, the form of the temperature dependence 
of <r sc , which reflects the topology of the SC gap, changes as a function of pressure. The most 
striking change is in the low-temperature behaviour of cr sc (T). At ambient pressure cx sc (T) 
shows a flat behavior below T/T c ~ 0.4, whereas the high-pressure data exhibit a steeper 
(linear) temperature dependence of <r sc (T) below T/T c ~ 0.4. We show in the following how 
these behaviors indicate the appearance of nodes in the gap function. 

B. Temperature and pressure dependent magnetic penetration depth 

In order to investigate a possible change of the symmetry of the SC gap, we note that 
A(T) is related to the relaxation rate cr sc (T) by the equation [34] : 

= omm >^Ty « 

where 7 ^ is the gyromagnetic ratio of the muon, and $0 is the magnetic-flux quantum. 
Thus, the flat T-dependence of <r sc observed at p = 0 for low temperatures (see Fig. 2b) 
is consistent with a nodeless superconductor, in which A ~ 2 (T) reaches its zero-temperature 
value exponentially. On the other hand, the linear T-dependence of cx sc observed at p = 2.25 
GPa (see Fig. 2b) indicates that A ~ 2 (T) reaches A~ 2 (0) linearly which is characteristic of 
line nodes. This is the most striking result of this communication: Pressure in an optimally- 
doped Fe-HTS can tune a nodeless gap into a nodal gap. Although this qualitative analysis 
is robust, it does not elucidate whether these nodes arise due to a nodal s +_ state or a 
d -wave state. 

To proceed with a quantitative analysis, we consider the local (London) approximation 
(A £, where £ is the coherence length) and first employ the empirical a-model. The latter, 
widely used in previous investigations of the penetration depth of multi-band superconduc¬ 
tors [31 135W40] . assumes that the gaps occuring in different bands, besides a common T c , are 
independent of each other. Then, the superfluid density is calculated for each component 
separately [3] and added together with a weighting factor. For our purposes, a two-band 




model suffices, yielding: 


A- 2 U) A- 2 (r,A al ) , A- 2 (r,A„, 2 ) 

A -2 (0) 1 A -2 (0, A 0il ) + 2 A- 2 (0,A„, 2 )’ () 

where A(0) is the penetration depth at zero temperature, Ao i? ; is the value of the i-th SC 
gap {i — 1, 2) at T — 0 K, and cuj is the weighting factor which measures their relative 
contributions to A -2 (i.e. uq + u> 2 = !)• 

The results of this analysis are presented in Figs. 3a-f, where the temperature dependence 
of A~ 2 for Bao. 65 Rb 0 . 35 Fe 2 As 2 is plotted at various pressures. We consider two different 
possibilities for the gap functions: either a constant gap, Ao,; = A*, or an angle-dependent 
gap of the form Ao= Aj cos 2p>, where ip is the polar angle around the Fermi surface. 
The data at p = 0 GPa are described remarkably well by two constant gaps, Ai = 2.7(5) 
meV and A 2 = 8.4(3) meV. These values are in perfect agreement with our previous results 
[3] and also with ARPES experiments [29], pointing out that most Fe-based HTS’s exhibit 
two-gap behavior, characterized by one large gap with 2A 2 /kB,T c = 7(2) and one small gap 
with 2 Ai//cbP c = 2.5(1.5). In contrast to the case p = 0 GPa, for all applied pressures 
A ~ 2 (T) is better described by one constant gap and one angle-dependent gap, confirming 
the presence of gap nodes, as inferred from our qualitative analysis. Note that a fitting to 
two angle-dependent gaps is inconsistent with the data. 

To understand the implications of the fitting to a constant and an angle-dependent gap for 
finite pressures, we analyze the two different scenarios in which nodes can emerge: a nodal 
s +_ state (with gap functions of different signs in the hole and in the electron pockets) and 
a d- wave state. In the former, the position of the nodes are accidental, i.e. not enforced by 
symmetry, while in the latter the nodes are enforced by symmetry to be on the Brillouin 
zone diagonals. Schematic representations of both scenarios are shown in Fig. [4] where a 
density plot of the gap functions is superimposed to the typical Fermi surface of the iron 
pnictides, consisting of one or more hole pockets at the center of the Brillouin zone, and 
electron pockets at the border of the Brillouin zone. In this figure, we set the accidental 
nodes of the s + ~ state to be on the electron pockets, as observed by ARPES in the related 
compound BaFe 2 (Asi_ 3 ,P x )2 [17]. Note that in the d- wave state, while nodes appear in the 
hole pockets, the electron pockets have nearly uniform gaps. Thus, the fact that the fitting 
to the a-model gives a constant and an angle-dependent gap is consistent with a d- wave 
state. 
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Figure 4: Schematic representation of the nodal s” 1 and d- wave states. In both panels, a 
density plot of the gap function is superimposed to a representative Fermi surface consisting of a 
hole pocket (h) at the center and an electron pocket (e) at the borders of the Brillouin zone. In the 
nodal states (upper panel), the nodes are not enforced by symmetry (here they are located at 
the electron pockets). In the d -wave state (lower panel), the nodes are enforced by symmetry to 
be on the diagonals of the Brillouin zone, and therefore can only cross the hole pockets. 
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To contrast the scenarios of a nodal s +_ gap and a d- wave gap, we consider a microscopic 
model that goes beyond the simplifications of independent gap functions of the cc-model 
discussed above. In this microscopic model, the fully coupled non-linear gap equations 
are solved for a hole pocket h and two electron pockets e\^, and the penetration depth is 
calculated at all temperatures. The free parameters are then the density of states of the 
pockets, the amplitude of the pairing interaction, and the gap functions themselves (details 
in the supplementary material). For simplicity, the anisotropies of the electron pockets are 
neglected, the Fermi velocities of the pockets are assumed to be nearly the same, and the gaps 
are expanded in their leading harmonics. Thus, for the nodal s +_ state we have A^ = A^ 0 
and A ej = A ei0 (r ± cos2<^ e ), whereas for the d- wave state it follows that A h = A/^o cos 2ph 
and A ei = ±A ei o- Note the difference in the position of the nodes in each case: while for 
the d,- wave case they are always at ph = ±7t/4, for the nodal s" 1 the nodes exist only when 
r < 1 at arbitrary positions p e = ±|arccos r. The results of the fittings for the pressures 
p = 1.57 GPa and p =2.25 GPa imposing a nodal state are shown in Figs. 3c and 
f. Remarkably, we find in both cases that the best fit gives r —> 0. This extreme case 
is, within our model, indistinguishable from the fitting to the d- wave state, since in both 
cases the nodes are at p = ±7t/4 (albeit in different Fermi pockets). We note that from 
the fits one cannot completely rule out the possibility of small but non-vanishing values of 
r. Therefore, at least within our model, a nodal state is compatible with the data only 
if the accidental nodes are fine-tuned to he either at or very close to the diagonals of the 
electron pockets for a broad pressure range. Since the position of the accidental nodes is 
expected to be sensitive to the topology of the Fermi surface, and consequently to pressure, 
it seems more plausible that the gap state is d- wave, since in that case the position of the 
gaps is enforced by symmetry to be along the diagonals of the hole pockets regardless of the 
value of the pressure. 

The pressure dependence of all the parameters extracted from the data analysis within 
the a model are plotted in Figs. 5(a-c). From Fig. 5a a substantial decrease of A (0) with 
pressure is evident. At the maximum applied pressure of p = 2.25 GPa the reduction of 
A (0) is approximately 15 % compared to the value at p = 0 GPa. Both Ai and A 2 show 
a small reduction upon increasing the pressure from p = 0 to 1.17 GPa, while above p = 
1.17 GPa the gaps values stay constant. On the other hand, the relative contribution 
of the small gap to the superfluid density increases by approximately factor of 2 for the 
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Figure 5: (Color online) Pressure dependence of various quantities of Bao. 65 R-t>o. 35 Fe 2 As 2 . 
(a) The magnetic penetration depth A(0) and A -2 (0). (b) The zero-tempearture gap values Ai^O). 
The relative weight ui 2 of the small gap to the superfluid density. The dashed lines are guides to 
the eyes, and the solid lines represent linear fits to the data. 


maximum applied pressure of p — 2.25 GPa (see Fig. 5c), indicating a spectral weight shift 
to the smaller gap. The parameters extracted from the microscopic model are discussed in 
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the supplementary material. 


III. DISCUSSION 

The most essential finding of our paper is the observation that pressure promotes a nodal 
SC gap in Ba 0 . 65 Rbo. 35 Fe 2 As 2 . This conclusion is robust and model-independent, as it relies 
on the qualitative change in the low-temperature behavior of AA -2 from exponential to 
linear in T upon applied pressure. To our knowledge this is the first direct experimental 
demonstration of a pressure induced change in the superconducting gap structure in a Fe- 
ffTS’s. Two possible gap structures could be realized at finite pressures: a nodal s +_ state 
and a d -wave state. In the first case, the change from nodeless s +_ to nodal s +_ is a crossover 
rather than a phase transition SU023, whereas in the latter it is an actual phase transition 
that could harbor exotic pairing states, such as s + id j2TH23] or s + d [ 33] . 

Additional results provide important clues of how pressure may induce either an s H or 
a d -wave state. In the closely related optimally-doped compound Bao. 6 K 0 . 4 Fe 2 As 2 , Raman 
spectroscopy [26], as well as theoretical calculations [201 EU, reveal a sub-dominant d -wave 
state close in energy to the dominant s +_ state. Pressure may affect this intricate balance, 
and tip the balance in favor of the d -wave state. On the other hand, theoretical calculations 
have shown that the pnictogen height is an important factor in determining the structure of 
the s +_ SC order parameter [18]. A systematic comparison of the quasiparticle excitations 
in the 1111, 122, and 111 families of Fe-HTS’s showed that the nodal s +_ state is favored 
when the pnictogen height decreases below a threshold value of ~ 1.33 A [33], Hydrostatic 
pressure may indeed shorten the pnictogen height and consequently modify the s + ~ gap 
structure from nodeless to nodal. Although our fitting of the penetration depth data to 
both a microscopic model and an effective a-model suggest that the d -wave state is more 
likely to be realized than the nodal s +_ state, further quantitative calculations of the pressure 
effect are desirable to completely discard a nodal s + ~ state. 

Besides the appearance of nodes with pressure, another interesting observation is the 
reduction of A(0) under pressure, despite the fact that T c remains nearly unchanged. Inter¬ 
estingly, in the compound BaFe 2 As 2 _ ;r P a ;, a sharp enhancement of A(0) is observed as optimal 
doping is approached from the overdoped side ra. which has been interpreted in terms of a 
putative quantum critical point (QCP) inside the SC dome [351 - I37] . In Bao. 65 Rbo. 35 Fe 2 As 2 , 
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if such a putative QCP is also present, pressure is likely to move the system away from the 
putative QCP, which, according to the results of BaFe 2 As 2 _ a ;P x , would explain the observed 
suppression of the penetration depth at T = 0. This scenario does not explain why T c stays 
nearly constant under pressure, but this could be due to the intrinsic flatness of T c around op¬ 
timal doping in Bai_ a ,Rb a; Fe 2 As 2 . Note that a similar behavior for A (0) and T c with pressure 
has been recently observed in LaFeAsOi_ a .F x [IB], but interpreted in terms of the interplay 
between impurity scattering and pressure. To distinguish unambiguously between these two 
scenarios, pressure-dependent studies of the quasiparticle mass in Ba 0 . 65 Rbo. 35 Fe 2 As 2 are 
desirable, in order to probe whether a putative QCP is present or not in this compound. 


IV. CONCLUSIONS 

In conclusion, the zero-temperature magnetic penetration depth A (0) and the temper¬ 
ature dependence of A -2 were studied in optimally doped Bao. 6 sRbo. 35 Fe 2 As 2 by means of 
/i SR experiments as a function of pressure up to p ~ 2.25 GPa. The SC transition temper¬ 
ature stays nearly constant under pressure, whereas a strong reduction of A (0) is observed, 
possibly related to the presence of a putative quantum critical point. Our main result is the 
demonstration that in the investigated Fe-based superconductor a nodal SC gap is promoted 
by hydrostatic pressure. Model calculations favor a d -wave over a nodal s +- -wave pairing as 
the origin for the nodal gap. The present results offer important benchmarks for the elucida¬ 
tion of the complex microscopic mechanism responsible for the observed non-universaltiy of 
the SC gap structure and of high-temperature superconductivity in the Fe-HTS’s in general. 


V. METHODS 

Sample: Polycrystalline samples of Bao. 65 Rbo. 35 Fe 2 As 2 were prepared in evacuated 
quartz ampoules by a solid state reaction method. Fe 2 As, BaAs, and RbAs were obtained 
by reacting high purity As (99.999 %), Fe (99.9%), Ba (99.9%), and Rb (99.95%) at 800 
°C, 650 °C and 500 °C, respectively. Using stoichiometric amounts of BaAs or RbAs and 
Fe 2 As, the terminal compounds BaFe 2 As 2 and RbFe 2 As 2 were synthesized at 950 °C and 
650 °C, respectively. Finally, samples of Bai_ ;r Rb a ;Fe 2 As 2 with x = 0.35 were prepared from 
appropriate amounts of single-phase BaFe 2 As 2 and RbFe 2 As 2 . The components were mixed, 
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pressed into pellets, placed into alumina crucibles, and annealed for 100 hours under vacuum 
at 650 °C with one intermittent grinding. Powder X-ray diffraction analysis revealed that 
the synthesized samples are single phase materials. 

Pressure cell: Pressures up to 2.4 GPa were generated in a double wall piston-cylinder 
type of cell made of MP35N material, especially designed to perform //SR experiments under 
pressure [491150] . As a pressure transmitting medium Daphne oil was used. The pressure was 
measured by tracking the SC transition of a very small indium plate by AC susceptibility. 
The filling factor of the pressure cell was maximized. The fraction of the muons stopping in 
the sample was approximately 40 %. 

/iSR experiment: The measurements were performed using high-pressure //SR tech¬ 
nique, where an intense high-energy = 100 MeV/c) beam of muons is implanted in the 
sample through the pressure cell. In a //SR experiment nearly 100 % spin-polarized muons 
/i + are implanted into the sample one at a time. The positively charged // + thermalize at 
interstitial lattice sites, where they act as magnetic microprobes. In a magnetic material 
the muon spin precesses in the local field B M at the muon site with the Larmor frequency 
Up = 7 At /(27r)R At (muon gyromagnetic ratio 7 M /(27r) = 135.5 MHz T _1 ). By means of //SR 
important length scale of superconductor can be measured, namely the magnetic penetra¬ 
tion depth A. When a type II superconductor is cooled below T c in an applied magnetic 
held ranging between the lower (H&) and the upper (H c2 ) critical held, a vortex lattice is 
formed which in general is incommensurate with the crystal lattice, and the vortex cores will 
be separated by much larger distances than those of the unit cell. Because the implanted 
muons stop at given crystallographic sites, they will randomly probe the held distribution of 
the vortex lattice. Such measurements need to be performed in a held applied perpendicular 
to the initial muon spin polarization (so called TF configuration). 

Analysis of TF-//SR data: Our zero-held //SR experiments (see supplemental material) 
reveal a magnetic fraction of about 10 % in the sample, caused by the presence of diluted 
Fe moments as discussed in previous //SR studies. The signal from the magnetically ordered 
parts vanishes within the hrst 0.2 /is. Thus, the hts of TF data were restricted to times t > 
0.2 //s for all temperatures. 
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The TF p SR data were analyzed by using the following functional form: 


P(t) = A s exp 


i° 2 sc + rtmW 


+A pc exp 




cos (7 p B int)S t + <p) 
cos (7 p B intjPC t + ip). 


( 3 ) 


and A pc denote the initial assymmetries of the sample and the pressure cell, respectively. 
7/(27 t) ~ 135.5 MHz/T is the muon gyromagnetic ratio, tp is the initial phase of the muon- 
spin ensemble, and B int represents the internal magnetic field at the muon site. The relax¬ 
ation rates <r sc and cx nm characterize the damping due to the formation of the vortex lattice 
in the SC state and of the nuclear magnetic dipolar contribution, respectively. In the anal¬ 
ysis cr nm was assumed to be constant over the entire temperature range and was fixed to 
the value obtained above T c , where only nuclear magnetic moments contribute to the muon 
relaxation rate cr. The Gaussian relaxation rate cr pc reflects the depolarization due to the 
nuclear magnetism of the pressure cell. It can be seen from the FT’s shown in Figs, lc and 
d that the width of the pressure cell signal increases below T c . As shown previously |5T] . 
this is due to the influence of the diamagnetic moment of the SC sample on the pressure 
cell, leading to a temperature dependent <r pc below T c . In order to consider this influence, 
we assume a linear coupling between cr pc and the field shift of the internal magnetic field in 
the SC state: a pc (T) = cr pc (T > T c ) + C(T)(p 0 H intjNS - AD#mt,sc), where a pc (T > T c ) = 
0.35 ps -1 is the temperature independent Gaussian relaxation rate. poH int ^s and //o-f^int,sc 
are the internal magnetic fields measured in the normal and in the SC state, respectively. 
As indicated by the solid lines in Figs, la-d, the /iSR data are well described by Eq. (1). 
The solid lines in panels c and d are the FTs of the fitted curves shown in Figs, la and b. 
The model used describes the data rather well. 
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SUPPLEMENTAL MATERIAL 


A. Experimental Details 

1. Sample characterization 

The temperature dependence of the zero field-cooled (ZFC) and field-cooled (FC) dia¬ 
magnetic susceptibility of Bao. 65 Rbo. 35 Fe 2 As 2 measured in a magnetic field of fi 0 H = 1 mT is 
shown in Fig. 6(a). From the diamagnetic response the SC transition temperature T c is de- 



Figure 6: (Color online) (a) Temperature dependence of the zero-field cooled (ZFC )and field- 
cooled (FC) susceptibility % m obtained in an applied magnetic field of po H = 10 mT for op¬ 
timally doped Bao. 65 Rbo. 35 Fe 2 As 2 . (b) Specific heat C p /T as a function of temperature of 
Bao. 65 Rbo. 35 Fe 2 As 2 . The arrow denotes the superconducting transition temperature T c . 
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Figure 7: (Color online) (a) Specific heat jump A C p at the superconducting transition vs T c for 
Bao. 65 Rbo. 35 Fe 2 As 2 , plotted together with literature data for various FeAs-based superconductors. 
The line corresponds to AC P oc T 3 (after [1]). 

termined from the intercept of the linearly extrapolated zero-field cooled (ZFC) susceptibility 
curve with Xm = 0 line, and it is found to be T c = 36.8(5) K. The temperature-dependent 
heat capacity data for this sample plotted as C p /T vs T is shown in Fig. 6(b). The jump as¬ 
sociated with the SC transitions is clearly seen. Here the anomaly at the transition has been 
isolated from the phonon dominated background by subtracting a second order polynomial 
C Pjn fitted above T c and extrapolated to lower temperature. The quantity AC P /T with AC P 
= (C p - C p n ) is presented as a function of temperature in the inset of Fig. 6(b). Although 
there may be some uncertainty in using this procedure over an extended temperature range, 
the lack of appreciable thermal SC fluctuations, as evidenced by the mean-field-like form of 
the anomaly, means that there is very little uncertainty in the size of AC p . Bud’ko et. al. [I] 
found that in many ’122’ Fe-based superconductors the specific heat jump AC p at T c follows 
the empirical trend, the so-called BNC scalling AC p oc T 3 . This has been interpreted as ei¬ 
ther originating from quantum critically or from strong impurity pair breaking. A violation 
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of the BNC scaling was observed for Bai_ a; K. r Fe 2 As 2 for x > 0.7 [Tj and in addition a change 
of the SC gap symmetry was observed. The specific heat jump data for Bao. 65 Rbo. 35 Fe 2 As 2 
obtained in this work is added in Fig. 7 to the BNC plot taken from Ref. [Tj. Our data point 
lies perfectly on the BNC line. 

2. Muon spin rotation experiments 

Zero-field (ZF) and transverse-field (TF) pSR experiments at ambient and under various 
applied pressures were performed at the /iEl beamline of the Paul Scherrer Institute (PSI), 
Switzerland, using the dedicated GPD spectrometer. A gas-flow 4 He (base temperature ~ 4 
K) and a VARIOX cryostat (base temperature ~ 1.3 K) were used. High energy muons (pp 
= 100 MeV/c) were implanted in the sample. Forward and backward positron detectors 
with respect to the initial muon spin polarization were used for the measurements of the 
pSR asymmetry time spectrum A(t). The typical statistics for both forward and backward 
detectors were 6 millions. All ZF and TF pSR experiments were performed by stabilizing 
the temperature in prior to recording the /iSR-time spectra. Note that a precise calibration 
of the GPD results was carried out at the 7 tM 3 beamline using the low background general 
purpose instrument (GPS). The /iSR time spectra were analyzed using the free software 
package MUSRFIT [2]. 

In a /rSR experiment nearly 100 % spin-polarized muons /i + are implanted into the sample 
one at a time. The positively charged p + thermalize at interstitial lattice sites, where they 
act as magnetic microprobes. In a magnetic material the muon spin precesses either in the 
local or applied magnetic field B M at thepenetration depth A and the coherence length £. If 
a type II superconductor is cooled below T c in an applied magnetic field ranged between the 
lower (H c i) and the upper (H c2 ) critical fields, a vortex lattice is formed which in general 
is incommensurate with the crystal lattice and the vortex cores will be separated by much 
larger distances than those of the unit cell. Because the implanted muons stop at given 
crystallographic sites, they will randomly probe the held distribution of the vortex lattice. 
Such measurements need to be performed in a held applied perpendicular to the initial muon 
spin polarization (so called TF conhguration). 
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3. Results of the Zero-Field, pSR experiments 


It is well known that undoped BaFe 2 As 2 is not superconducting at ambient pressure and 
undergoes a spin-density wave (SDW) transition of the Fe-moments far above T c [3]. The 
SC state can be achieved either under pressure m or by appropriate charge carrier doping 
of the parent compound [6], leading to a suppression of the SDW state. Magnetism, if 
present in the samples, must be taken into account in the TF-/rSR data analysis. Therefore, 
we have carried out ZF-^SR experiments above and below T c to search for magnetism 
in Ba 0 . 65 Rbo. 35 Fe 2 As 2 . As an example, ZF-/iSR spectra recorded at T —5 K and 50 K of 
Bao. 65 Rbo. 35 Fe 2 As 2 are shown in Fig. 8a. There is no preccesion signal, indicating that 
there is no long-range magnetic order. On the other hand, we observed a significant drop 
of the asymmetry, taking place within 0.2 /rs. This is caused by the presence of diluted 
Fe moments as discussed in previous /iSR studies [7]. In order to quantify the magnetic 
fraction, the ZF-/rSR data were analyzed by the following function: 
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the first and the second terms describe the magnetic and nonmagnetic part of the signals, 
respectively. Aq is the initial asymmetry, fl is the magnetic volume fraction, and At (Al) is 



Figure 8: (Color online) (a) The ZF-^SR time spectra for Bao. 65 Rbo. 35 Fe 2 As 2 recorded above and 
below T c . The solid line represent the fits to the data by means of Eq. (1). (b) Temperature de¬ 
pendence of the magnetic fraction of Bao. 65 Rbo. 35 Fe 2 As 2 , extracted from the ZF-pSR experiments. 


26 









t (Lis) 


Figure 9: (Color online) ZF-//SR time spectra for Bao. 65 Rbo. 35 Fe 2 As 2 at various applied pressures 
recorded at the base temperature T = 1.4 K. The solid line represents the fit to the data by means 
of the sum of the Eq. (1) and a damped Kubo-Toyabe depolarization function to account for the 
pressure cell signal. 

the transverse (longitudinal) depolarization rate of the /iSR signal, arising from the magnetic 
part of the sample. The second term describing the paramagnetic part of the sample is the 
combination of a Lorentzian and a Gaussian Kubo-Toyabe depolarization functions M . a 
and A are the depolarization rates due to the nuclear dipole moments and randomly oriented 
diluted local electronic moments, respectively. The temperature dependence of the magnetic 
fraction obtained for Bao. 65 Rbo. 35 Fe 2 As 2 is plotted in Fig. 8b. The magnetic fraction at the 
base temperature was found to be only 8 %. Bearing in mind that the signal from the 
magnetically ordered parts vanishes within the first 0.2 ps in the whole temperature region, 
the analysis of transverse field data was restricted to times t > 0.2 ps. 

Figure 9 shows the ZF-^SR time spectra for Bao. 65 Rbo. 35 Fe 2 As 2 at various applied pres¬ 
sures. The ZF relaxation rate stays nearly unchanged between p = 0 GPa and 2.25 GPa, 
implying that there is no sign of pressure induced magnetism in this system. 
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B. Microscopic model for analyzing the penetration depth data of 


Bao.65R-bo.35Fe 2 As 2 


1. Model for pairing 


As a minimal model that accounts for the different superconducting states of the iron 
pnictides (nodeless s +_ , nodal s + ~, and d- wave), we consider a two-dimensional system with 
three isotropic Fermi pockets HD): one hole pocket h centered around V = (0, 0) and two 
electron pockets e\ and e 2 centered around = (tt, 0) and M 2 = (0,7r) (see Fig. 10). To 
describe the s + ~ state, the pairing interaction between the hole pocket h and the electron 
pocket e\ is assumed to be angular dependent with the form: 

V hei = V 0 (r - cos2(/>)hj(k)/r{(—k)e u (—p)eit(p) + h.c. , (5) 

where 0 is the polar angle measured relative to the center of the electron pocket, Vo is the 
interaction energy scale, and r is the relative amplitude of the angular-independent and the 
angular-dependent pairing interactions. Due to the tetragonal symmetry of the system, the 
pairing interaction between h and e 2 is: 

V he2 = Vo(r + cos20)h{(k)/r{(-k)e 2 ^(-p)e 2t (p) + h.c. . (6) 


Furthermore, to minimize the number of free parameters, we assume that the three 
pockets have the same Fermi velocity Vf, while the density of states can in principle be 
different pu/p e = P- Within this model, we obtain an s +_ state, where the SC gap of 
the hole pocket is a constant, A^, and the gap on the electron pockets is of the form 
A ei = A e (r — cos 20) and A e2 = A e (r + cos 20). Accidental nodes appear in the electron 
pockets if r < 1. Introducing the energy cutoff A c , we can write down the corresponding 
BCS-like gap equations: 
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where E ei ( k), E e2 ( k), and E h ( k) are the quasi-particle energy dispersions: 
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(8) 


E ei (k) = y/ef+ Af(r - cos20) 2 , E 2 { k) = \Je\ + A 2 e (r + cos20) 2 , E h ( k) 
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M 2 = (0, 7r) 



Figure 10: (Color online) Three pocket model used in our calculations. It is assumed that the 
system has one hole pocket h centered around T = (0, 0) and two electron pockets e\ and e 2 
centered around M\ = (vr, 0) and M 2 = (0,7r). 


To determine T c , we linearize the gap equations, yielding: 


de (3 (- 

A, = —A e p e V 0 (2r 2 + 1) / - tanh A- 

J o e ^ 

a a T r [ c de /3 c e 

A e = -A h 7jp e V 0 / — tanh — 

Jo e z 


PeYo ~ 


y/?/(2r 2 + 1) f — tanh —— 

Jo e 2 


-l 


To perform the fitting, we set T c to be fixed, and set the energy cutoff A c = 86meV (the 
results do not depend significantly on the choice of the cutoff). This provides a constraint 
on p e V o, rj, and r. When T < T c , the gaps are calculated based on the BCS Eqs. ([7]) and 


The expression for the penetration depth of a single-band system is: 



where / is the Fermi distribution function, e is the energy of the non-interacting system, 
and Ek is the quasi-particle energy dispersion. Applying this formula to our three pocket 
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model, we obtain 

,_ am v f['\,( 8 f 9/W ,/*• 

A mccp ^L Ae {aE h -arJ +p ‘ v 'L 


i [ d 0 2 

de / — cos i 
./ 2vr 


2 / A ' , f &<t> 2 

p e vj / de / —cos 

■/-A c 


2vr 


dE, 


A 2 (T) OCPefj 


2 + p 


(l - 2/(A c )) +Tj I de 


de e 

df_ 

dE h 


e2 

A, 


+ 2 


de 


df df 
dE ei de e 


d0 df 

2n dE P 


( 9 ) 


In the fittings, we will focus on the normalized penetration depth A 2 (T) /A 2 (0). 


2. Model for d-wave pairing 

To describe the d -wave superconducting state within our three band model, we consider 
the following form of the pairing interaction: 


V hei = V 0 (r - cos26»)h|(k)h{(-k)e u (-p)ei t (p) + h.c. 

V he2 = V 0 (r + cos 26 ) )h|(k)/r|(—k)e 2 4,(—p)e 2 t(p) + h.c. . 
where 9 is the angle around the hole pocket. The gap functions can then be written as: 

A ei = -A e2 = A e , A /l (k) = A h cos 29 . 


resulting in the BCS-likc gap equations: 


Ah — 2A e p e Vo 


'-A c 2 E e 

rA c 


de , $E e 
tanh- 


A e = A h r]p e V 0 / de 


'-Ac 


d9 cos 2 29 (3Eh 

-tanh- 

27t 2 E, 


h 


Here, r/ = p h /p e , E e = y/e | + A 2 e , and E h = \]e 2 + A 2 cos 2 29. Repeating the same steps 
as for the s +_ case, we obtain the penetration depth: 


A 2 (T) oc p e vj 


2 + r) 


(1 — 2/(A c )) + p 


de 


'0 


d 9 df 
2n dE h 


+ 2 


' A ' de 9/ 


'0 


dE P 


( 10 ) 


Comparing the expressions for the d -wave case to the expressions we derived for the s + 
case, Eqs. (j7|) and we note that they can be mapped onto each other if r = 0. In this 
extreme case, changing rjd — > d/r\ s , Vo,d —> pV o, s /2, and A h o A e leads to the same gap 
equations and penetration depth expression. With these replacements, both s and d pairing 
give the same A -2 (T)/A~ 2 (0). Therefore, we conclude that the penetration depth cannot 
distinguish between nodal-s +_ and d -wave if the nodal-s +_ is the extreme case with r = 0. 
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3. Fitting Results 


We now fit the experimental data A -2 (T) /A -2 (0) of optimally-doped Bai_ x Rb x Fe 2 As 2 
to find the values of p e V r 0 , r/, and r for different pressures. Note that the value of T c imposes 
another constraint on these three parameters, as explained above. Figs. 11a, b and c show 
the fitting for the s + ~ model for P = 0, P = 1.57 GPa, and P = 2.25 GPa, respectively. For 
the P = 0 case, we find equal gap amplitudes and no nodes, as seen by ARPES experiments 
in the related compound Bai_ x K I .Fe 2 As 2 . We see that the fitting is not as good in the region 
immediately below T c . We will discuss this issue in more details below. For the pressurized 
samples, the fitting is overall better and indicates a nodal state (r < 1). The value of 
the density of states ratio ph/Pe is little affected by pressure (as expected, since no charge 
carriers are introduced), and is consistent with the value of a nearly compensated metal. 

Surprisingly, the best fittings for both the P = 1.57 GPa and P = 2.25 GPa cases give 
r = 0, where the nodes on the electron pockets are fixed at 6 = ±7 t/4 . This is a very 
special case of the accidentally nodal s +_ state, since by symmetry there is no reason for r 
to vanish. To make this point more transparent, in Fig. 12a we plot the non-zero pressure 
data and the theoretical urves for the penetration depth for various values of r - keeping 
all the other parameters constant. Clearly, 0 < r < 1 gives worst fittings than r = 0. What 



Figure 11: (Color online) The temperature dependence of A~ 2 (T)/A _2 (0) measured at various 
applied hydrostatic pressures of Bao. 65 Rbo. 35 Fe 2 As 2 . The square symbols are experimental data 
and the red curves are the theoretical functions, (a) Fitting for the P = 0 data, which suggests a 
nodeless state, (b) and (c) Fitting for P = 1.57 GPa and P = 2.25 GPa. The fitting suggests that 
nodes exists on the two electron pockets at the angles 6 e = ±7r/4 and ±37r/4. 


31 








Figure 12: (Color online) (a) Effect of the electron pocket gap anisotropy on the penetration 
depth at p = 2.25 GPa. The electron gap is nodal if r < 1, and becomes nodeless if r > 1. The 
low temperature data clearly shows that the gap is nodal, but the data near T c seems to be better 
described by a nodeless state, (b) Fitting for the zero pressure case with the Fermi velocity ratio 
Vh/v e being a free parameter. The fitting improves with respect to Fig. 6a, but the values of Vh/v e 
and ph /Pe seem to be too large or too small. 

we also found is that r = 10 - i.e. a nodeless superconducting state - describes the data 
better near T c , on the expense of a very bad fitting at low temperatures - where the nodal 
behavior is evident. 

As we discussed in the previous section, a nodal-s +_ state with r = 0 is indistinguishable 

- for fitting purposes - from a d- wave state. Since there is no symmetry reason to have 
r = 0 in our simple model, or even r<l over a wide pressure range, we interpret this result 
as an indirect indication that a d- wave state is more likely to be the state of the pressurized 
samples. 

Finally, we comment on the difficulty of the fittings to capture the behavior near T c 

- particularly for the sample at ambient pressure (see Fig. 11a). One reason could be 
the presence of inhomogeneities, which would require a distribution of gaps to be taken into 
account, instead of a single gap value. Another reason could be related to our choice of fixing 
the Fermi velocities to be the same for both the electron and hole pockets. To investigate this 
possibility, we lift this restriction and allow Vh/v e to also be a fitting parameter. The result 
is shown in Fig. 12b. Clearly, we obtain a better fitting, but not only p e V r 0 is relatively large, 
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but the ratios phjp e and Vh/v e are very large or very small, which is difficult to reconcile with 
the Fermi surface of these materials. Most likely, additional pockets are necessary to capture 
the full temperature dependence of the penetration depth. Nevertheless, our microscopic 
model provides results that agree with those obtained from the a-model fitting, particularly 
in the low-temperature regime, suggesting that a d -wave state is more likely to be realized 
than a nodal s +_ state. 
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